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Abstract 

The zero temperature phase diagram of XXZ spin chain in external magnetic field 
is investigated at low energies using path integral approach.lt has been shown by spin 
wave analysis and then by nonlinear sigma model transformation that below some crit- 
ical field the system undergoes quantum Kosterlitz-Thouless phase transition. Above 
that critical field for low anisotropics the system is saturated while at strong anisotropy 
the system has Ising antiferromagnetic order. 

1 Introduction 

Considerable effort has been made to understand the long range (low energy) 
behaviour of quantum spin systems specially in low dimensions. The impor- 
tance of these studies is in their applicability in various physical phenomena 
in for example quantum Hall ferromagnets[l], high superconductors [2] and 
many other magnetic systems. On the other hand we know very few models 
to be exactly solvable, nevertheless development of field theoretical approaches 
has greatly changed the language in this part of strongly correlated systems 
and has provided us with effective methods to be applied to our problems of 
spin systems. One of the great observations by means of field theoretical ap- 
proach has been the problem of ground state of ID Heisenberg antiferromag- 
net: after separating out the rapidly varying fluctuations we obtain a nonlin- 
ear sigma model(NL(TM) with a topological term which is responsible for de- 
stroying the long range order of ground state for integer spin chains according 
to Haldane's conjecture[3] while we know that this problem is exactly solv- 
able by Bethe's ansatz only for spin half problem [4]. The Haldane's conjecture 
has been now confirmed by strong experimental and numerical evidence [4, 5, 6]. 
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The method of separating out short wavelength fluctuations which resulted in 
NLfjM+topological term for antiferromegnetic spin chains is now widely used 
for various spin systems including Fcrrimagncts[7] and spin laddcrs[8].In this 
paper I am going to apply this method of obtaining low energy behaviour to 
anisotropic spin chains. 

2 Spin Wave Analysis 

For the anisotropic Heisenberg model in magnetic field: 

1 ^ 

^ = ^ E [2 (^'^^^ + ^^^/) + ^^''^^] + E B • s (1) 

<ij> j = l 

where < A < 1 there are numerical evidence [9] that under some critical value 
for the magnetic field Be the Sz fluctuations become negligible and the sys- 
tem behaves as XY model in long range. On the other hand wc know that 
the isotropic Heisenberg antiferromagnet becomes massive by turning on mag- 
netic field [8]. These two considerations suggest that Be <x (1 — A)"" for XXZ 
model. This may be examined by some crude spin wave analysis: I Apply the 
Holstein-Primakov transformation which is some kind of perturbation around 

saturation point: Sf = s — aja^ and S~ = (2s)^a|(l — -^)^ ■ Here cre- 
ates a boson to lower the 5^ by one.We are in a situation where the num- 
ber of bosons are low enough to approximate the square root and obtain: 
S~ ~ (2s) 2 a^. Applying this transformation to the Hamiltonian (1) for B = Bz 
and Fourier expanding we find: 

n^J2l^ijk-zX)Js-B]alak (2) 

k 

in which 7^ = '^^^e^'''^ — 1-k- a represents the translation vectors to the 
neighboring sites and z is the number of nearest neighbors. 
In one dimension the spin wave spectrum becomes gappful: Ek = 4Js(cosafc — 
X) — B with the gap A = 4 Js(l — A) — B.This estimation confirms my suggestion 
and sets u = 1. Under Be the saturation point becomes unstable. This elementary 
analysis will be invalid as B approaches zero and is valid just near Bg.This result 
shows that near Be the correlation length of Sz fluctuations £^ (x\ B — Be ■ 
In section 3 after separating out short range Neel field we can see that the 
correlation length of Sz fluctuations is actually ^oc(l — A)~2 forA<l. 
For the case of A > 1 the analysis is a bit complicated. In this regime to find the 
fluctuations around Ising antiferromagnetic order(IAF) we can again apply the 
Holstein-Primakov transformation for two sublattices: Sf = s — ajaj and S^ = 

\/2sa|(l — "27^ )^ for i in sublattice A and for sublattice B: Sf = —s + blbt and 

5r = ^/2i(l-^)i6i. Here we have two kind of Bosons one a^b^) lowers(adds) 
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the Sz by one in sublattice A(B).By applying the first order approximation: 
Ri \/2saJ for A and ~ \/2sbi for B, and inserting in the Hamiltonian (1) 
we find: 

= Js^jk{akb-k + alb^-k) +" E ^fc^*^ + ^ E ' 

k k k 

Where a = JsXz + B and (3 = JsXz — B and we have Fouric;r c'xpandcd the 
operators in the last equality . Now to diagonalize (3)1 apply the (canonical) 
Bogoliubov transformation : 

Cfc = Mfeflfe - Wfe&Lfc (4) 

dk = Ukbk - Vka^_f, 

For the transformation to be canonical we must have:u^ — = 1 and for (3) to 
be diagonal we must have: Js^{v? + u^) + (a + (3)uv = 0. No problem arises if 
we assume u and v are real. Finally we obtain: 

w«E4'^4c/fe+E4''4^fe (5) 

k k 

4^'^^ = 2Js{\^ - cos^ afc)i ± B. 

One of The branches is always gapfuU which has to be the spectrum of fluctu- 
ations trying to align spins against the magnetic field direction and the other 
one is also gapfuU which tries to saturate the system, but in this branch the gap 
vanishes at: Be = 2Js{X^ ~ 1)^ • This result shows that for A > 1 the correlation 
length ^ oc (A — 1)~2 that is v = ^ above the critical point A = 1. This result 
is more reliable for the following reason: 

we can compute the site magnetization reduction due to quantum fluctuations 
in our spin wave approximation: 

|A<5,>|= vldk= /°( .^"^^ -bdk (6) 

a a VA2 - 1 ■ 

Where F is the complete elliptic integral of first kindWe can see in Fig. I that 
this reduction is finite at A > 1 and vanishes very fast as A grows. 

Now we can sketch the line of transition from saturation to lAF and other 
phases as follows(fig.II): 

Bc(A) = 2JsVA2 - 1 A > 1 (7) 
B^(A) = 4Js(l - A) A < 1. 
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3 Field Theoretic Discussion 



In this section I am going to analyze the anisotropic spin chain by means of 
nonlinear sigma model transformation. The Minkowskian action for anisotropic 
Heisenberg chain W = J2f=i ~ SjLi B • S with the symmetric rank 

two tensor a is : 

Sm = s'S^Swzinj]- dxoy2o:^„n':'n'■_^_l - dxoV'B-n. (8) 

Where xq is the real time and /? = l/fcsT that in this paper we discuss the zero 

temperature situation jSh — > cx^.The first term known as Wess-Zumino term is 
a Berry phase for each spin and represents the area of a cap of the unit sphere 
formed by the closed path passed by the unit vector n(a;o) from a;o = to 
Xo = /37i[10] .Under the critical magnetic field Be we can separate out rapidly 
varying fluctuations as follows: 

First of all we stagger the field: {—lYiij then we write rij as a short 

range Neel field Ij plus fluctuations around it mj : nj = uij +ao(— 1)^1^ where 
= 1, m • 1 = and ao is the lattice constant. This method works only for 
two sublattice systems close to Neel state which is not the eigenstate of H of 
course! 

Expanding the action up to second order in 1 and integrating it out we obtain 
Z = J Dm5(m2 - 1) e'^'/f where : 

Seff= / (fx[—aaum''m'' + '^a^^m'^dlm'"] + ^ [ cfxm- {doui x aim) + (9) 
J ao 2 2 J 

' f - 1 
- / a~^[(m X doni)u,{in x doni)^](fx - — - — B • (m x dom). 

J JsXao 
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For the XXZ model a^^ = Js^ diag{l, 1, A) by which we find: 
^eff = ^(dom.f - ^{dim.f - £(1 - A)m^ - ^[{drm.f + (10) 
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-2{mymzdomzdomy + mxruzdornxdoniz + A ^mxniydomydomx)} + 

B 9 
— — (m X dom) ■ z + — e„,.m • {du.m x di,va) . 

JsAao oTT 

where for the last term n^v = 0,1. The last term is the Pontryagin index (or 
winding number) of field m as a map from 5^ (real space-time) to S"^ (con- 
figuration space) and is a topological term and 9 = 27rs[10]. This is the well 
known 0-term that for isotropic Heisenberg chain resulted in Haldane's conjec- 
ture, however in XXZ chain as we can see in (10) due to massive fluctuations 
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of niz (A oc 1 — A) the z component will be negligible in ground state, in other 
words this component is short range and so at long wavelengths(low energies)is 
irrelevant. To find the low energy behavior we can use the mean-field approxi- 
mation rUz « 0. In this approximation the winding number will be ineffective. 



4 Topological Phase Transition at T=0 

In our mean-field approximation removing the z component in the Lagrangian 
density results in: 

>Ce// = -^[{dim^r + idimyf] + ^Ij^lidomyfml + (11) 
+{domx)'^my - 2mxmydomxdomy] - — {rrixdomy - mydorrix) 

Where in the last line we have parametrized m as (cos (p, sin (p, 0) and scaled the 
space-time so that the spin wave velocity is one and g = (2A)2 is the coupling 
constant. This is the Lagrangian for massless fluctuations of XXZ model in 1-1-1 
dimensions. (The term coupled to magnetic field is a total derivative, i.e. the 
magnetic field has no effect on in plane excitations in mean-field approxima- 
tion.) The Lagrangian in (11) has topological defects for the field (p{xo,x) is 
an angle. These defects drastically alter the behavior of system when they are 
present. In this scenario(Kosterlitz-Thouless), at weak couplings {g < gc) the 
ground state has bound vortices (defects in space-time) and is partially ordered 
with algebraically falling correlations. At strong couplings the vortices become 
unbound and the partial order completely vanishes and the correlations fall off 
exponentially with correlation length ^ oc 1/ ln(25)).This drastic change is a 
second order quantum phase transition at zero temperature without symmetry 
breaking. 



5 Summary 

The low energy phase diagram of antiferromagnctic XXZ chain in magnetic 
field is shown in fig. II . At < A < 1 and under Be ~ 4J,s(l — A) we have 
two kind of ground states: at A < \kt the ground state is partially ordered 
while at A > \kt the ground state has free vortices and is disordered. Above 
the critical magnetic field the z component is no longer ignored and we have 
saturated ferromagnetic(SF) order. At A > 1 and B < Be obviously the system 
has antiferromagnctic Ising (lAF) behaviour at low energies while in this regime 
and above Be we have again saturated phase. For B = it is inferred that we 
have a quantum phase transition at A = 1 from disordered phase to ordered 
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lAF phase which has aheady been shown in continum hmit for spin half case 
by Bosonization method[10]. 
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Figure caption 

Fig. I) Site magnetization reduction due to quantum fluctuations at ground state 
in the regime of higher anisotropy A > 1 . 

Fig. II)The low energy phase diagram of the XXZ antiferromagnetic chain. The 
straight line shows the threshold of saturation. The curve ^/B'^ + 1 is the thresh- 
old of Ising Antiferromagnetism(IAF). 
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